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SUMMARY 


Measurements are presented for pressure gradients induced by a 
laminar boundary layer on a flat plate in air at a Mach number of 9.6 
and for the drag of thin wings at a Mach number of about 6.8 and zero 
angle of attack. The pressure measurements at a Mach number of 9.6 were 
made in the presence of substantial heat transfer from the boundary 
layer to the plate surface. The measured pressure distribution on the 
surface of the plate was predicted with good accuracy by a modification 
to insulated-plate displacement theory which allows for the effect of 
the heat transfer and temperature gradient along the surface on the 
boundary-layer displacement thickness. 

The total drag of thin wings with square and delta plan forms was 
measured at a nominal Mach number of 6.8 over a reasonably wide range 
of Reynolds numbers. The total drag was found to be greater than can 
be explained by adding a classical value of laminar skin friction to the 
estimated pressure drag. The difference is, in general, explained by 
the increase in skin friction (20 to 40 percent) caused by the boundary— 
layer-induced pressures . 


INTRODUCTION 


The considerable distortion of the flow field about plane surfaces 
or slender bodies in hypersonic flow due to boundary-layer displacement 
effects is now well known qualitatively. This phenomenon is a result of 
the low mass involved in laminar boundary-layer flow and can achieve 
considerable importance at high Mach numbers and low Reynolds numbers. 

The most readily observable result of the presence of the boundary layer 
is a local increase in surface pressure which appears in conjunction 
with an increase in skin friction and heat transfer. A considerable 
body of quantitative theoretical evaluations of this effect is now avail- 
able. These solutions are concerned mainly with the sharp-leading-edge 
problem and can be considered to fall into two classes, the so-called 
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"weak interaction" and "strong interaction" solutions. Expressed rather 
simply, the weak-interaction regime is one in which the self-induced 
pressure gradient does not have an important effect on the boundary-layer 
growth, whereas in the strong-interaction regime large induced pressures 
occur which do have an important effect on the boundary-layer growth and 
necessitate consideration of the mutual interaction between the boundary 
layer and the boundary-layer-induced pressure gradients. 

The parameter which basically determines the region in which inter- 
action effects can be important was first given by Lees and Probstein 
(ref. l) . This parameter takes into account Mach number, Reynolds num- 
ber, and temperature level of the external stream. Other evaluations of 
the weak- and strong-interaction regions have been made by Bertram, Shen, 
Li and Nagamatsu, Stewartson, Lees, Kuo, and Pai and Shen (refs. 2 to 9) • 

The experimental work in air has been mainly supplied by Bertram 
for a Mach number of 6.9 (refs. 2 and 10) and Kendall for a Mach number 
of 5.8 (ref. 11). The main drawback to this data taken in air has been 
that only the lower end of the range in which the self-induced effects 
become important has been explored. The data which have been obtained 
in helium by Bogdonoff and Hammitt at the Princeton helium hypersonic 
tunnel (refs. 12 and 13 ) cover a range in which the self-induced effects 
can be of major importance. Most of these helium data, however, include 
a considerable leading-edge-bluntness effect which makes difficult the 
comparison with viscous-interaction theory. 

Among the boundary-layer-induced effects are the pressure gradient 
induced by the presence of the boundary layer and the increase in skin 
friction associated with this induced pressure gradient. With the 
recent availability of a nominal Mach number 10 nozzle for the Langley 
11 -inch hypersonic tunnel there was the possibility of considerably 
extending the range of the interaction data in air. This nozzle has 
been utilized to obtain new measurements of the boundary-layer-induced 
pressures, although heat transfer from the boundary layer to the model 
is present in these experimental results. These measurements are pre- 
sented in this report, together with measurements of the total drag at 
zero lift of wings with thin square and delta plan forms in the nominal 
Mach number 7 nozzle of the 11-inch hypersonic tunnel. These total-drag 
measurements are utilized to obtain an indication of the increased skin 
friction that can be ascribed to the induced pressure field. 


SYMBOLS 


A 


constant in Blasius formula for skin friction 


a 


fraction of root chord of delta wing 
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a',b' skin-friction constants in viscous-interaction theory 

B second virial coefficient in equation of state for imperfect 

gas 

b wing semi span 

C - M-y-^oo / m-oo^w 

C,D constants in power formulas for turbulent skin friction in 

appendix D (eqs. (D12) and ( Dl4) ) 

C 1 = ^l/^l^w 

wave- or pressure-drag coefficient at zero angle of attack 
Cd^t total drag coefficient at zero angle of attack 

C D t leading-edge drag coefficient based on wing area 

Op average skin-friction drag coefficient 

Cp average skin-friction drag coefficient without displacement 

effect 

Cf local skin-friction drag coefficient 

Cf local skin-friction drag coefficient without displacement 

effect 

C-n mew pitot pressure coefficient 

JU y iILcLX 

Ch + average pressure coefficient over leading edge of wing 

c local chord 

c effective chord of a wing for the determination of average 

skin friction (see appendix D, eq. (Dl)) 

Ac length of sting contained within the wing plan form (fig. 3) 

c r root chord 

d sting thickness (fig. 3) 
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F = 0.149(7 - 1) (7 + 1) 

f length of sting base flare (fig. 5) 


g length of sting outside of wing plan form (fig. 5) 

2 


L plate length 

M Mach number 


m,n,q 


Ni 


Kn 


reciprocals of exponents in power laws for skin friction in 
appendix D 


Knudsen number, 


Molecular mean free path 
Leading-edge thickness 


Np r Prandtl number 


P 

Pc 

Po' 

R 

R 



Rrp 


R t 


R x 

S 

T 


T 


c 


pressure 

critical pressure of gas 

total pressure measured by pitot tube 

Reynolds number based on undisturbed free-stream conditions 
gas constant 

Reynolds number with root chord of delta wing as character- 
istic length 

Reynolds number for transition based on distance from leading 
edge in free-stream direction 

Reynolds number with leading-edge thickness as characteristic 
length 

Reynolds number with x as characteristic length 
local lateral dimension of plate or wing from center line 
absolute temperature 
critical temperature of gas 
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Tp oo static temperature for constant value of 7 

t leading-edge thickness 

u velocity 

X,Y,Z nozzle coordinates (fig. l) 

x longitudinal distance from leading edge of plate 

y distance measured laterally on model from center line 

a angle of attack 

(3 Sutherland constant 

7 ratio of specific heat at constant pressure to that at 

constant volume 

6 boundary -layer thickness 

x 

6 boundary -layer displacement thickness 

e semiapex angle of delta wings 

p dynamic viscosity 

p gas density 

X viscous-interaction parameter 

X c viscous-interaction parameter with Reynolds number based on 

r root chord of delta wing 

a) exponent in power law for viscosity 

Subscripts : 

e recovery conditions 

ins insulated wall 

0 supply conditions 

turb turbulent flow 

w 


wall 
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1 edge of boundary layer 

A delta wing 

oo free stream 

Superscript: 

* ratio of conditions taking variation in y into account to 

conditions with constant y 

APPARATUS AND METHODS 


Tunnel and Nozzles 

The investigation was conducted in the Langley 11-inch hypersonic 
tunnel, which operates intermittently. Three different nozzles were 
used during the tests: 2 two-dimensional nozzles, each of which provided 

a Mach number of slightly less than 7* and a three-dimensional nozzle 
with which a Mach number close to 10 could be obtained. The first two- 
dimensional nozzle had contours machined from steel and was replaced 
after the test program had started by a nozzle having contours constructed 
qP invar . Invar was used for the contour plates of the second nozzle 
because thermal gradients in the nozzle blocks caused objectionable 
deflection of the first minimum of the steel nozzle. In addition, the. 
Invar nozzle was designed so that pressure gradients normal to the hori- 
zontal plane of symmetry were minimized. The three-dimensional nozzle 
had a square test section and first minimum, though nowhere. else was 
the cross section square. This design allows the installation of. con- 
ventional schlieren windows and reduces the possibility of incurring the 
axisymmetric imperfections at the wall that might occur in a circular 
nozzle and which would tend to focus along the axis of the nozzle. 

A description of the tunnel may be found in reference 14 and a 
description of the steel nozzle and a calibration at a stagnation pres- 
sure of 25 atmospheres (R per inch = 250,000) in reference 15.. Addi- 
tional calibration data for the steel nozzle and some calibration data 
for the Invar nozzle are contained in reference l6. Figure 1 presents 
some of the calibration data obtained in the three-dimensional nozzle. 


Models and Supports 

Of the wings investigated, two had square plan forms (fig. 2) and 
the rest were a family of delta wings with symmetrical double-wedge 
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sections, a thickness ratio of 2 ~ percent, and semiapex angles varying 

from 30° to 8 °. The surfaces were ground and the leading edges were 
from 0.001 to 0.002 inch thick. The designations of these wings and 
their dimensions and sting dimensions are shown in figure 3 - 

For the measurement of surface pressures the model shown in fig- 
ure 4 was used. This model was also equipped with thermocouples for 
the determination of skin temperature. The leading-edge thickness of 
this model was close to 0.0005 inch. 


Ins trumen tati on 

The drag of the wing models was measured with a two-component strain- 
gage balance. This balance is temperature compensated and its sensi- 
tivity to uneven heating has been reduced to tolerable limits by shielding 
and insulation. For a more complete description of this balance see 
reference 17 , where it is designated as the "sensitive two-component 
force balance." 

The base pressure for use with the sting corrections to the wing 
force measurements and the surface pressures for the model shown in 
figure 4 were measured by means of the aneroid-type six-cell recording 
unit described in reference 14. The stagnation pressure was measured 
with Bourdon gages with an accuracy of about one-half of 1 percent. 


Surface-Film Flow Studies 

V 

Surface flow studies of the square -plan- form wings and wings 3* 

4-B, and 4-C were made by photographing the patterns made by streaming 
graphite and fluorescing mineral oil under ultraviolet light during a 
run. The wings were coated with SAE 30 lubricating oil before the run 
and graphite was spotted along the leading edge. The cameras were 
equipped with suitable filters to photograph the fluorescing oil to best 
advantage . 


Tunnel Conditions 

For the force tests of the wings at a Mach number near 7 the tunnel 
was operated at a stagnation temperature near 1 , 130 ° R, and for the 
square -plan-form wings it was operated through a stagnation-pressure 
range from 7 to 39 atmospheres (with Reynolds number per inch in the 

range from 0.08 x 10° to 0.43 X 10°). Calibrations of these nozzles 
show a slight effect of Reynolds number per inch in the test range. 
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For example, at Reynolds numbers per inch of 0.08, 0.14, 0.20, 0.26, 

and 0.40 (x 10^) the corresponding Mach numbers are 6.67, 6-77; 6.82, 

6.86, and 6.90. 

For the determination of pressure distributions in the three- 
dimensional nozzle the stagnation temperature was close to 1,600° R and 
the flat surface of the pressure model (fig. 4) was located parallel to 
the XY-plane of the nozzle at Z = 3/4 inch, with the leading edge at 
X = -2 inches (fig. l) • The stagnation pressure was varied through the 
limited range from 31 to 46 atmospheres (with Reynolds number per inch 
6 6 

from 0.075 x 10° to 0.11 x 10° when the variation in 7 is taken into 
account). The Mach number in this nozzle (fig. l) is about 9-6 if caloric 
imperfections (variable 7) are taken into account (ref. 18) . 


PRECISION OF DATA 


For the force data obtained in the Mach number 7 nozzles, random 
errors in coefficients arise mainly in reading and recording balance 
forces and stagnation pressures, whereas the source of systematic errors 
is mainly in the possibility of an incorrect assessment of Mach number 
and the effect of the sting pressure field on the wing. The possible 
error resulting from errors in Mach number determination is small in 
comparison with the possible random errors, being about one-half of 
1 percent of the coefficient values at Reynolds numbers per inch greater 
than about 200,000 and rising to perhaps 0.8 of 1 percent below this 
value . 

The maximum possible errors in drag coefficient (ACp) due to randan 
factors, including sting corrections, are as follows (Reynolds number is 
based on root chord for the wings with delta plan form) : 


Wing 

Plan form 

R 

ac d 

5-percent-thick 
diamond sections 

Square 

0.3 X 10 6 

•5 

•9 

1.4 

0.00075 

.00044 

.00024 

.00016 

Wedge - s lab-wedge 
sections 

Square 

• 5 

• 5 

• 9 

1.4 

2.4 

.00031 
. 00031 
.00021 
.00015 
.00014 
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Wing 

Plan form 

R 

ac d 



0.7 X 10^ 

0.00051 

1 

Delta 

1.4 

.00035 



2.1 

.00027 



*9 

.00047 

2 

Delta 

1.4 

. 00029 



2.1 

.00025 



1.0 

. 00051 

3 

Delta 

1.6 

.00031 



2.7 

.00024 



1*5 

.00046 

4 

Delta 

2.0 

.00036 



2.7 

. 00032 



2.3 

.00023 



2.6 

.00020 

4-A, 4-B, 4-C 

Delta 

3.0 

4.0 

.00017 

.00014 



4.8 

.00009 



5.4 

.00014 


The pressures measured on the flat plate at M - 9*6 have a possible 
error of about 2 percent due to inherent instrument inaccuracies. 


RESULTS AMD DISCUSSION 


Theoretical Considerations and Comparison With Previous Data 

The considerable distortion of the flow field about* plane surfaces 
or slender bodies in hypersonic flows due to boundary- layer displacement 
effects is now well known qualitatively. In addition, a considerable 
body of quantitative theoretical evaluations of this effect is now avail- 
able. These solutions in general fall into two classes: the so-called 

"weak interaction" and "strong interaction" solutions. Expressed rather 
simply, the weak-interaction solution is one in which the self-induced 
pressure gradient is assumed not to affect the boundary -layer growth, 
whereas in the strong-interaction solution account is taken of the effect 
on the boundary-layer growth of the mutual interaction between the 
boundary layer and the boundary-layer-induced pressure gradients. The 


10 


NACA TN 4133 


parameter which basically determines the region in which interaction 
effects can be important was first given by Lees and Probstein (ref. 1) . 
This parameter is commonly designated X, where 


X = 




As shown subsequently, the weak-interaction theory for the insulated 
case applies in the region of X from 0 to approximately 1, whereas 
strong-interaction theory applies in the region of X » 1. A leading 
edge of infinitesimal thickness is postulated. 

Insulated plate in weak-interaction regime . - Into the category of 
weak-interaction solutions fall the theories formulated by Lees and _ 
Probstein (ref. 1) and Bertram (ref. 2) . In terms of the parameter X 
the pressure distribution and skin friction for the insulated flat plate 
are as follows for large values of M and small flow deflections. 

The Lees-Probstein solution for pressure distribution, to the 
second order, is 

-2- = 1 + 0.5997(7 - 1)X + 0.0907(7 - 1) 2 (7 + 1)X 2 (Np r = l) (1) 

-E_ = 1 + 0.5667(7 - 1)X + 0.077 7(7 - 1) 2 (7 + 1)X 2 (N pr = 0.725) (2) 

1*00 


If Bertram's expression for the shape of the boundary-layer displacement 
thickness is put into the power-series expansion for pressure (with the 

product small), to correspond with the Lees-Probstein theory, 

the second-order result is 


■£- = 1 + 0.6387(7 - 1 )X + 0.1027(7 - 1) 2 (7 + DX 2 (n P i . = l) ( 3 ) 


The difference in coefficients between equations (l) and (3) arises from 
a slight difference in the boundary-layer thickness values used in the 
theories. First-order weak-interaction theory consists of unity plus 
the first term in X in equations (l) to (3). 
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A simple expression for the local skin-friction coefficient is 
obtained in reference 2 from flat-plate theory, for a linear velocity 
profile in the boundary layer with no heat transfer: 


C f = 


2^0 


u n 


Pcn U < 


2 6 


(4) 


An expression for 6 can be obtained by taking the result from refer- 
ence 2 to hypersonic approximation (after Von KArmAn, ref. 19)* If the 
variables in the resulting equation are evaluated at the local conditions 
at the edge of the boundary layer, the result is: 

2.555 L -~- Mp /Cp x 

6 = 2 (5) 


Substituting equation ( 5 ) into equation (4) and utilizing the perfect gas 
law and the adiabatic relations between pressure, density, and tempera- 
ture gives for local skin-friction coefficient the expression 


c f & = 0.7854^/^ (6) 

v c 

or, for comparison with the flat-plate theory without pressure gradient, 

= = \/p/p» (7) 

c f 

In a fashion similar to that in which equation (6) was obtained, equa- 
tion (5) may be written in terms of undisturbed free-stream conditions, 
and the normalized boundary-layer thickness is 


6^R 

x Vc 


= 2.555 


M 


, 2 <p/pJ - 1/2 


2 


( 8 ) 
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As will be seen by comparison with the results of strong- inter act ion 
theory presented subsequently, equations ( 6 ) to ( 8 ) will also apply as 
an approximation through the strong- interaction regime. 

Kuo's closed- form results in reference 8 apparently apply to the 
weak-interaction region of plates since, with M -»», his equation for 
pressure distribution is almost identical to equation (l) . Kuo's evalu- 
ation of local skin-friction coefficient is composed of two opposing 
terms . With M -» <» the result obtained is 


Although equation ( 9 ) indicates an initial increase in skin friction 
Over the no-interaction case with increasing X, the skin friction 
eventually reaches a peak and then decreases rapidly. 

Insulated plate in strong-interaction regime .- The strong-interaction 
regime has been investigated by Lees (refs. 6 and 7 ), Bertram (ref. 2), 
and Li and Nagamatsu (ref. 4) . This regime is characterized by a boundary- 
layer growth of the form 



where 



00 


and 


F = 0.149(7 - 1) (7 + 1) 
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for large values of M and X. (See, for example, p. 26l of ref. 6.) 
The classical flat-plate laminar -boundary- layer growth used in the 
weak-interaction solution to obtain the induced-pressure distribution 
is 


where M is large. 

The results of the strong-interaction solutions apply, in general, 
for values of the interaction parameter X much greater than 1. Two 
main approaches to the strong-interaction problem are available. One 
is Shen's assumption (ref. 3) > later extended by Li and Nagamatsu 
(ref. 1+) , that the edge of the boundary layer is also the shock induced 
by the boundary layer, and the second is the approach by Lees (refs. 6 
and 7 ) £in d Bertram (ref. 2) in which the edge of the boundary layer (or 
the displacement thickness) is taken as the boundary of a new body. The 
effect of these differences in concept are discussed by Lees in refer- 


The results of the Li-Nagamatsu theory may be put into the following 
form (for an insulated flat plate), which is consistent with the hyper- 
sonic approximations used in the theory: 


ence 6. 



C f^ = 2 27 (7 - D 3A *1/2 


- 7 = ^-^1 77 

^ (7 + l) lA 



*2 > 3-91 

(7 + l) (7 - 1) 
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= = 1 + ( 4.12 - O.506 - 1 )X 

c F V V7 - 1/ 


(^X 2 < 3-91 \ 

\ (7 + l) (7 - l) / 


For the Lees strong-interaction theory applied to the insulated 
flat plate the following results are obtained: 


_P_ 

Pc* 


0.52X + 0.92 


(First order; 7 = 7/5; Npr = l) 


JL = 0.92X (Zeroeth order; 7 = 5/3; Np r = l) 

Poo 


hi!, 0.55 * 1/2 

{c 

0.83X 1 / 2 

Cf 

^ - 2.20X 1 / 2 
1.66X 1 / 2 

Cp 


(Zeroeth order; 7 = 7/5; N Pr = l) 


In Bertram's theory, effects approaching strong interaction can be 
taken into account by an iteration procedure described in reference 2 . 
However, with this method there is no distinction between strong and weak 
interaction, since the iteration procedure gives a continuous variation 
between the two regimes. Pai and Shen (ref. 9) present a theory in 
which, by matching the solution to both the strong- and weak-interaction 
regimes, a result is obtained which is also applicable through both 
regimes . 

Shown in figure 5 are the pressure distributions given by the various 
theories for an insulated plate with 7 = 7/5> together with Kendall's 
experimental results (ref. ll) for a flat plate in thermal equilibrium 
in air (surface temperature essentially equal to recovery temperature). 
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From various considerations , including the Knudsen number and the dis- 
tance of the first orifice from the leading edge (in terms of leading- 
edge thickness) , the effect of leading-edge blunting on these pressures 
is believed to be small. Leading-edge-thickness effects are discussed 
further in appendix A. The curve for Bertram's theory was calculated 
for M = 9.6 by the iteration scheme described in reference 2 by use 
of isentropic compressions and continuous variations in dM/dx. Though 
this Mach number is considerably higher than 5-8, the curve should be 
applicable to the data for lower Mach numbers since this curve agreed 
with the results from the iteration process for M = 6.86 presented in 
reference 10 in the region of X where they overlapped. The process 
of using the result of a first iteration to start a second iteration, 
and so on, does not appear to be convergent. However, by averaging 
between assumed and resulting pressure distributions for successive 
iterations the results can be made convergent in the range of X below 
about 10. The second-order weak-interaction_theory for air appears to 
agree well with experiment in the range of X from about 0 to 2. The 
Lees strong-interaction, the Pai-Shen, and the Bertram theory agree well 
with experiment over the entire range in which experimental data are 
available (up to X « 3) . The Li-Nagamatsu theory would not be expected 
to apply at the low values of X. The curve shown in figure 5 for the 
Li-Nagamatsu theory does not agree with that given in figure 17 of refer- 
ence 12. However, it can be shown that the curve given in reference 12 
is incorrect and does not apply even when C = 1. 

The theoretical local and average skin-friction coefficients for 
the insulated plate sire shown in figure 6 divided by the value given 
by flat-plate theory without interaction effects. Also shown are the 
experimental average skin-friction coefficients over the limited range 
of X obtained by Kendall (ref. ll) , using the Eimer-Nagamatsu results 
as the interaction-free denominator (ref. 20). There is essential 
agreement between the Bertram theory with pressure gradient and the l£es 
strong- interact! on theory; however, the Li-Nagamatsu theory predicts 
skin frictions which are considerably higher than those of either of 
the foregoing theories. The experimental results apparently favor the 
Bertram theory, which at the largest value of X of the experimental 
data predicts a skin friction about 25 percent below that given by the 
Li-Nagamatsu theory. Such correlation between theory and experiment, 
however, cannot be considered to be conclusive in this case since some 
question arises as to the accuracy of the data. If the local skin- 
friction coefficients are multiplied by the factor 2, which is assumed 
to give the average skin-friction coefficient, Cp^R 83 1.07 from the 
Eimer-Nagamatsu experiments whereas theoretically, for the no- interaction 
case, C f /R is approximately 1.35 (ref. 21). The cause for the dis- 
crepancy is unexplained. In addition, although the investigations of 
references 11 and 20 were conducted in the same wind tunnel, the methods 
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of evaluating skin friction were different. This difference is expected 
to he an additional source of error, though less important than the dis- 
crepancy involved above. 

N oninsula ted plate .- The case of the noninsulated plate has not 
been investigated as extensively as that of the insulated plate. Quali- 
tatively, it is easily seen that heating the plate above the recovery 
temperature will increase the boundary-layer-interaction effects through 
thickening of the boundary layer, whereas the reverse will be true if 
the plate is cooled below recovery temperature because the boundary layer 
will be thinned. 

In treatments of the interaction problem when heat transfer is 
present, the wall temperature is generally considered constant. For the 
weak-interaction case, Lees and Probstein give a solution (ref. 1) 
including heat transfer from which the induced-pressure effect may be 
obtained in closed form. For the strong-interaction region there are 
several recent analyses of the effect of heat transfer when the wall 
temperature is invariant. Lees (ref. 6) gives zeroeth-order results 
when the heat transfer is very large (wall temperature near 0°) and Li 
and Nagamatsu (ref. 22 ) give zeroeth-order results for a wide range of 
wall temperatures and 7 values of 1.4 and 1.67. These results agree 
in regions where the two may be compared, even though reference 22 makes 
the considerably simplifying assumption that C = 1 in the linear formula 
for viscosity. However, the values given in reference 22 for the coeffi- 
cient used to evaluate the Stanton number are too high by a factor of 2 . 
The results of Pai and Shen (ref. 9), as mentioned previously, apply 
through the weak- and strong- interaction regimes. 


Experiment 

Pressure measurements in air at = 9.6 . - The pressure distribu- 

tions at a Mach number of 9.6 were obtained on a flat plate in the three- 
dimensional nozzle with heat transfer from the boundary layer to the 
plate surface. The heat transfer is a result of the short duration of 
the test run combined with the initially low flat-plate temperature 
(about room temperature) and the high stagnation temperature (about 
1 600° R) . This pressure distribution is shown in figure 7, together 
with the values given by theory for an insulated flat plate (from 
fiK. 3) • The measured pressures are considerably overestimated by 
insulated-plate theory. The pressures presented are those for 85 seconds 
after the start of a run and the temperatures for the same time are shown 

in figure 8. 

As a first approximation it was decided to calculate the effect of 
the temperature gradient and heat-transfer conditions on boundary-layer 
thickness independently of the pressure-gradient effects. To this end 
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the Chapman-Rubesin theory (ref. 23) was utilized with a free-stream 
Mach number of 9*6. The second-degree equation shown at the top of 
figure 8 was taken to represent the temperature distribution (the 
recovery temperature was arbitrarily selected to be about O.85 of the 
stagnation temperature) . Computations were then carried out to determine 
the displacement thickness at the 15-, 30- , 60-, and 100-percent-chord 
points. The results of these calculations are shown in the two lower 
plots of figure 8. The corrections to the insulated-plate results were 
obtained simply by modifying the slope of the boundary-layer displace- 
ment thickness calculated for the insulated flat plate, as follows: 

d 6 * 

d6 x _ 6 X d-Sjns + ^jns 5jns 

~ " ^ ne ^ ai 

The ratios 5 X /b x ns and d(b x /5 x ns ) / d ( x / L ) were obtained from the 

Chapman-Rubesin theory, and the values for d6^ n s / ^ ^ ins/ ^ were 

computed according to the theory with pressure gradient of reference 2. 

This correction to the slope of the boundary- layer displacement 
thickness on an insulated flat plate resulted in the pressure distribu- 
tion shown in figure 7 as the dashed line. The agreement with the 
experimental results is surprisingly good. For comparison with the 
computations that include heat transfer and temperature gradient, the 
effect of a constant ratio of wall to stream temperature was computed 
from the Chapman-Rubesin theory by using an approximately average tem- 
perature ratio from the experiments (T w /Too « 9) • The resulting value of 
6 X /B?n S is shown in figure 8; it is not an average value of the change 

in the displacement thickness due to heat transfer but, as expected, is 
too low a value and thus would underestimate the pressures on the average. 
The heat transfer itself, however, has much more effect on the pressures 
than the temperature -gradient term. 

In figure 7 the data at the lowest values of (highest 

Reynolds numbers) approach the base line of zero pressure difference at 
a faster rate than is predicted by theory. This is apparently the same 
phenomenon as was observed in the Princeton helium tunnel on similar 
models (figs. 11, 12, and 14 of ref. 12) and eliminated in the final 
plots of the helium- tunnel results (figs. 15 to l8 of ref. 12). There 
can apparently be an under or over pressure, depending on the mounting 
and the shape of the back and side of the model, and figure 11 of refer- 
ence 12 indicates that there is a strong Reynolds number effect. These 
side and base effects appear to be propagated for large distances both 
upstream and laterally. 
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In the present tests pressure orifices were located laterally at 
the middle and last station on the plate (see fig. 4) . The resulting 
pressure differences are shown in figure 9 &s a function of the local 
lateral distance to the orifice nondimens ionali zed with respect to 
local lateral distance to the model edge. There is considerable scatter 
in the data, much of which is attributable to the low pressure being 
measured (about 1 mm of mercury) but the general trend of a decrease in 
pressure as the model edge is approached is evident. The decrease in 
pressure from the center to the outermost orifice is considerably larger 
for the rearmost set of orifices. 

All the data from orifices on or near the model center line have 
been included in figure 7; thus some of the data are from the orifices 
between 2 and 4 inches from the leading edge (X less than about 2.5). 
This is the region in which the side effects discussed previously are 
important, and therefore the results from this region may be somewhat 
in error. 

Figure 10 presents the experimental results originally given in 
figure 7; however, instead of modifying the theory to account for heat 
transfer and temperature gradient, the theoretical increment in pressure 
attributable to heat transfer has been added to the data for M = 9-6. 
The good agreement of these data with the l£es, Pai-Shen, and Bertram 
theories is notable, as is the correlation with Kendall's data for 
M = 5.8. 

Real-gas effects .- In the evaluation of Mach number and the other 
associated flow properties of a gas it is important to determine the 
effects associated with actual gas flows. In most supersonic wind 
tunnels it is sufficient to use the perfect gas law and assume that y 
is constant. At the high supply pressures and stagnation temperatures 
at which hypersonic tunnels are operated these assumptions, however, 
can be invalidated. 

An analysis of the effect of variable specific heats on the gas 
properties has been given by Eggers (ref. 24) . The equations given in 
reference 24 have been combined with the Sutherland formula to determine 
the effect of caloric imperfections on parameters basic to the boundary- 
layer analysis. The additional assumption has been made that the static 
temperature is much less than the stagnation temperature. For the range 
of stagnation temperature shown in figure 11 (1,000° R to 5,000 R) , 
this assumption corresponds to the restriction that the static tempera- 
ture be less than approximately 500° R- The equations which describe 
the curves in figure 11 are given in appendix B. For the present tests 
the stagnation temperature was close to 1,600° R. For a pressure ratio 

of 30 X 10" 6 (see fig. l) the Mach number was found to be close to 9-6. 
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The Reynolds numbers computed for the conditions at this Mach number 
were corrected according to figure 11 with (3/Tp^ = 2.1. In the present 
tests the effects of variable gamma were not large. The question then 
arises as to the magnitude of imperfect-gas effects. 

With the assumption that the Berthelot equation of state in virial 
form is applicable, equations describing imperfect-gas effects have been 
derived in appendix B. In equation (Bll) the following values were then 
used: T c = 238.5° R, p c = 547 lb/sq in. abs (the values for air); 

T q = 1,600° R, P 0 = 647 lb/sq in. abs, p^ = 0.0194 lb/sq in. abs 
(experimental values); T OT = 82° R, 7 = 7/5 (assumed values). The 
resulting static-temperature ratio given by equation (Bll) was found 
to be only one-half of 1 percent greater than would be given by the 
ideal-gas equation (p = pRT, 7 = 7/5) . This increase is small compared 
with the 3.3-percent increase in this same temperature parameter indi- 
cated by equation (B3) for the effect of variable 7. Therefore imperfect- 
gas effects were ignored in analyzing the data. 

Drag measurements at ^ = 6.8.- The results of total-drag measure- 
ments in coefficient form for two square wings and a series of delta 
wings are presented in figure 12. Leading-edge thickness was determined 
for each wing from an examination of the leading edge through a calibrated 
microscope. The pressure acting on the leading edge was assumed to be 
given by the average between C r ^ may and the average pressure coefficient 

over the front half of a cylinder. The drag coefficient given by the 
average pressure coefficient (Cp,t = (5/6)Cp,max) was subtracted from 

the total drag coefficient for presentation in figure 12. For the delta 
wings C! was assumed to vary as the square of the sine of the 

P y IIIcLa. 

semiapex angle. Following are the estimated leading-edge drag coef- 
ficients in terms of projected wing area: 


Plan form 

Wing 

t, in. 

C D,t 


Diamond sections 

0.001 

O.OOO38 

Square 

Wedge - s lab -wedge 
sections 

.002 

.00056 


1 

0.001 

0.00018 


2 

.001 

.00006 

Delta 

3 

.001 

.00002 


4 

.002 

.00002 


4-A, 4-B, 4-C 

.0008 

.00000 
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For the purpose at hand, that of determining skin friction by using 
wings with a very low pressure drag, the square wing with 5-percent- thick 
diamond sections had an undesirably high pressure drag; however, this 
wing was the only one for which local pressures were available (the 
investigation reported in ref. 17) • In the final analysis the square 
wing with wedge-slab-wedge sections was hampered by a drag due to leading- 
edge thickness which was about double the value that would be obtained for 
inviscid pressure drag if the leading-edge thickness were nonexistent 
(Cp,p = 0.00030 compared with values given in preceding table). 

For the square wings (at the top of fig. 12) the drag obtained by 
adding the skin friction, with the effect of the self-induced pressure 
gradient included, to the theoretical inviscid pressure drag is generally 
closer to the experimentally determined drag than if the classical flat- 
plate skin friction is used. The total experimental drag is 10 to 
30 percent higher than is given by a combination of classical laminar 
skin friction and the inviscid pressure-drag coefficient. There are a 
great many unknowns still involved in the effects of viscosity on the 
drag of the rear surface, including tip, trailing-edge, and displacement 
effects. According to pressure measurements on the center line of the 
square wing with diamond sections at zero angle of attack (fig. 4(a) of 
ref. 17) the displacement effect on the front surfaces causes an incre- 
mental increase in pressure drag of 11 percent; however, this increase 
is partially compensated by an incremental decrease of about 40 percent 
in pressure drag of the rear surfaces due to displacement and trailing- 
edge shock effects. The result is an indicated 28-percent decrease in 
pressure drag from the theoretical prediction. 

An indication of tip effects on the square wing with the diamond 
section, at least at one Reynolds number, can be seen in figure 13(a) 
from an oil-flow study of the wing surface at a Reynolds number of 

1.2 x 10^. The oil flow on this wing was examined by means of both 
still and motion pictures. The significant extent of tip effects appears 
to be confined to the rear surface and consists of a roughly triangular 
region indicative of either a shock or separated flow. The dark areas 
near the leading edge of the wings shown in figures 13 and 14 are caused 
by the high laminar shear near the leading edge, which quickly wipes the 
oil from this region. 

Figure 13(b) shows the oil flow on the upper surface of the square 
wing with wedge -slab-wedge sections, and figure 13(c) shows the flat 
lower surface of this wing. The rear surface of this wing is small, 
and although the flow phenomenon is similar to that on the diamond- 
section wing, the separation is confined to a correspondingly smaller 
portion of the trailing wing tip. In addition to this separation (or 
shock) a three-dimensional edge flow can be seen. On the lower surface 
the region of three-dimensional edge flow is delineated by the high shear 
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which scrubbed the oil from the surface near the side edges. On this 
surface the angle that this scrubbed area makes with the wing edge is 
about 2°. On the upper surface a disturbance can be seen to extend 
farther inboard than the area actually scrubbed clean. The total area 
affected seems to make an angle with the wing side edge of about 3 « 5 °* 

For both wings the area significantly influenced by the tip flow is 
small; however, the effects on the total drag may be significant because 
of the indicated increase in surface shear in the tip region. 

At the bottom of figure 12 are the data for the family of 2^ - percent- 

thick delta wings with semiapex angles varying from 30 ° to 8 ° and sym- 
metrical double-wedge sections. The theoretical skin-friction-coefficient 
curves for these wings differ from those presented for the square wings 
only in their adaption to the delta plan form, as shown in appendix C. 

The wave drag shown in the figure for these delta wings is that given by 
two-dimensional shock-expansion theory for the streamwise wing section 
(assuming the compensating pressure effects on the front and rear wing 
surfaces discussed previously for the square wings) . This wave drag has 
been added to the skin friction with displacement effect to obtain the 
curve coded with short dashes in this figure. For the range of semiapex 
angles shown there is good correlation of • the experimental data over the 
entire range of Reynolds number (based on root chord) . At the lower 

Reynolds numbers ^R Cp < 1-3 X 10^ the theoretical skin friction with 

displacement effect plus the two-dimensional wave drag gives an excellent 
representation of the data (wings 1 , 2 , and 3 ) and is 20 to 30 percent 
higher than the drag obtained by a combination of the classical skin 
friction and inviscid pressure drag. Above this Reynolds number, however, 
there is a progressive departure of the experimental data from the laminar 

theory. This result for R n > 1.3 X 10^ is taken to be indicative of 

u r 

transition of the boundary layer. The data for the square wings in fig- 
ure 12 indicate a similar trend at the higher Reynolds numbers. It 

should be noted that the value of 1.3 X 10^ for the Reynolds number based 
on root chord corresponds to a different unit Reynolds number for each 
wing. For delta wings 1, 2, and 3 this Reynolds number corresponds to 

Reynolds numbers per inch of 0.30 X 10^, 0.24 X 10^, and 0.19 X 10^, 
respectively. 

Oil-flow studies of the delta wing (fig. 14) have shown that the 
only deviations from disturbance-free flow on the surface are manifested 
in the region of the sting. The disturbance is indicated by the dark, 
approximately triangular area emanating from the point of the sting, 
where the oil has been wiped from the surface. Wing 1 (e = 30°) shows 
very little sting interference. (See fig. 14(a).) This wing and 
wings 2 and 3 had a very narrow sting (0.091 to 0.092 inch) . In con- 
trast, the oil-flow photographs of wings 4-C and 4-B (figs. 14(b) 
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and 14(c)) show well-defined regions of sting interference. The wings 
of the number 4 series had thicker stings (see fig. 3) than wings 1, 2, 
and 3 . The results for wings 4-A, 4-B, and 4-C in figure 12 indicate 
that the pressure field from the stings has only a small direct influence 
on the drag, even though the length of sting covering the root chord was 
varied from 0.43 of the roof chord down to 0.19 of the root chord. In 
view of the low pressure drag of these airfoils, this result is not 
surprising. For wing 4-C, which was tested in two nozzles with very 
nearly the same test-section Mach number, there was no really significant 
difference between the data obtained in the two nozzles. The data from 
the two nozzles is, in general, within the band given by the possible 
error. (See section entitled "Precision of Data.") 

Computations were made to determine the effect that the sting would 
have if it were assumed to cause transition of the boundary layer. To 
simplify matters, a triangular flat plate with zero inclination to the 
stream was assumed, with transition parallel to the leading edge at 
various fractions of the root chord, a. For comparison with the hypo- 
thetical wing, the experimental results from figure 12 for the delta 
wings are shown in figure 15 with an assumed constant value of pressure 
drag coefficient (0.00037) subtracted. The curve for transition fixed 
at about 0-7 of the root chord fits the data over a considerable range 
of Reynolds numbers. Because of the relatively small increase in 
measured drag over the drag predicted when classical flat-plate laminar- 
boundary-layer theory is used, attributing the increase to displacement 
effects must be considered inconclusive. However, from figure 12 it 
appears that this increase in drag, for the lower Reynolds numbers at 
least (considering the thinness of the mounting sting and the thickness 
of the boundary layer), is caused by displacement effect. At Reynolds 

numbers greater than about 1.3 X 10^ the drag is higher than can be 
accounted for by any cause but transition. 


CONCLUDING REMARKS 


Measurements are presented for boundary-layer-induced pressure 
gradients on a flat plate in air at a Mach number of 9.6 and for the 
drag of thin wings at a Mach number of about 6.8 and zero angle of 
attack. The pressure measurements at a Mach number of 9-6 were made 
in the presence of substantial heat transfer from the boundary layer 
to the plate surface. The measured pressure distribution on the sur- 
face of the plate was predicted with good accuracy by a modification 
to insulated-plate displacement theory which allows for the effect of 
the heat transfer and temperature gradient along the surface on the 
boundary-layer displacement thickness. 
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The total drag of thin square and delta wings was measured at a 
nominal Mach number of 6.8 over a reasonably wide range of Reynolds 
numbers. This total drag was found to be greater, by as much as 30 per- 
cent, than can be explained by adding a classical value of laminar skin 
friction to the estimated pressure drag. The difference is, in general, 
explained by the effect of boundary-layer-induced pressures on skin 

friction. At Reynolds numbers above about 1.3 X 10^ an additional 
increase in skin friction was measured, which is attributed to transi- 
tion of the boundary layer. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 15, 1957* 
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APPENDIX A 


LEADING-EDGE-THICKNESS EFFECTS 


The question of the contribution of leading-edge effects to the 
measured pressures arises. Two obvious factors which have to be eval- 
uated are the (characteristic) boundary-layer thickness in comparison 
with the leading-edge thickness and the mean free path of the gas in 
comparison with the leading-edge thickness. 

The assumption is made that this boundary- layer thickness for the 
laminar case may be represented by the simple equation 


The Reynolds number with radius of curvature of the edge of the boundary 
layer as the characteristic length is then 


At the leading edge 


With a linear velocity profile and M -»» on an insulated plate. 


6 = K 


x 


(Al) 





where C is the coefficient in the viscosity relation 


r 

4G 
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or, to eliminate the viscosity ratio as such, the power law or Sutherland ! s 
law for viscosity may be used to give, with hypersonic approximation for 
the insulated plate and Np r = 1, 


C = 




(o)-l) 





+ 


p 


where a) is the exponent in the power law for viscosity. Thus 

1-Kjo 

/y - l\ 2 1+0) 

K = 2.553G-— — ) 


and 


R t,s ■ 6 -5 18 ( 21 T 1 ) 


2(1-Kjq) 


or 




^2(1 w 


R+ 


(A2) 


The question of whether the flowing gas will "see” the leading 
edge fully, imperfectly, or not at all should be answered by the 
Knudsen number Nj^, which is the ratio of the molecular mean free 

path to the characteristic length (in this case the leading-edge 
thickness) : 

„ _ R L _ 1 . 256^7 M 

“ R t " R t 

where Rp is the Reynolds number with mean free path as the character- 
istic length and the slip-flow regime is taken as 


loo < < 
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The fact that equation (Al) does not represent the boundary-layer 
growth near the leading edge is ignored because for present purposes all 
that is desired is an index of the boundary-layer-induced effects on the 
plate. Whether equation (A2) can give some idea of the relative magni- 
tude of the boundary-layer-induced effects and leading-edge-thickness 
effects can perhaps be decided by recourse to the available reports of 
experimental work. The work of Bertram (ref. 10) , Kendall (ref. ll) , 
and Hammitt and Bogdonoff (ref. 13) is useful in this regard. The range 
of these tests is shown in the table below: 


Author 

Ref. 

Ho 

R t,5 

R t 

R t,6 

R t 

Rt 

NlCn = Rt 

Bertram 

10 

6.9 

625 

370 to 1,960 

1.7 to O.32 

0.028 to 0.0052 

Kendall 

11 

5.8 

260 

14 to 47 

19 to 5-5 

0.62 to 0.18 

Hammitt and 
Bogdonoff 

13 


11.8 

11.4 

3,660 

3,300 

J940 to 42,500 

3.9 to 0.07 

0.02 to 0.0004 

Present 

tests 

— 

9.6 

900 

38 to 56 

24 to 16 

O.38 to 0.26 


For the experiments of Kendall (ref. ll) , the results of which are 
probably the best yet obtained in air, the combination of the large 
values of the ratio g/R-t shown in the table and a large Knudsen 
number which is apparently well into the slip-flow regime indicates that 
the leading-edge-thickness effects should be small compared with the 
boundary-layer effects. These data did not show an effect of R^ on 
surface pressure in the range of R t of the investigation. In contrast, 
the range of R^ 5/Rt for Bertram's tests in air (ref. 10) indicates 

that the leading-edge-thickness effects should range from small to a 
magnitude approximately equal to that of the boundary -layer effects and 
the leading-edge thickness is in general in the continuum regime. The 
results of pressure measurements in reference 10 show leading-edge - 
thickness effects over the test range that are of the order expected 
from consideration of equation (A 2) of the present paper. 

Hammitt and Bogdonoff 's test range in the Princeton helium tunnel 
(ref. 13) is indicated (in the preceding table) to extend from a region 
in which the boundary-layer-induced effects are about 3 times the leading- 
edge-thickness effects to a region where the boundary-layer-induced 
effects are small compared with the leading-edge -thickness effects; how- 
ever, over the entire range the leading edge is in continuum flow. The 
pressure data of reference 13 bear out this prediction in that, for values 
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of the index s/^t less than approximately 0 . 2 , the pressure data 

correlate well when plotted against the inviscid parameter x/t. Sig- 
nificant deviation from correlation is first noted at R^g/R^. = 0.55> 

and at R^g/R^ = 0.9 the measured pressure rises are about twice the 

values given by the essentially inviscid pressure data at the same value 
of x/t. The deviation from correlation is even greater at 

*t,6/Rt = 3.9. 

The data of the present tests are in the same region of slip flow 
for the leading-edge thickness as Kendall's data, and by equation (A2) 
the viscous effects are still larger in comparison with leading-edge- 
thickness effects. The present results have values of the ratio R^ 5 /^t 

which are as large as the largest value encountered in the Princeton 
tests (ref. 13 ) and the leading edge is farther into the slip-flow 
region. Leading-edge effects would not be expected in the present 
experiments and none could be detected over the limited range of R^. 
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APPENDIX B 


EFFECT OF VARIABLE 7 AND GASEOUS IMPERFECTIONS 


The problem of evaluating the effect of caloric imperfections 
(variable specific heats) may be greatly simplified by considering the 
stagnation temperature to be much larger than the static temperature . 
In this case the analysis of reference 24 can be used in the following 
simplified form (perfect gas law) . 


For the Mach number, 



(Bl) 


where 0 is the molecular vibrational constant . 


With the values of p 
required temperature ratio is given by 


and Pq obtained experimentally, the 


"0 



0 e 


; 9/T 0 


*/*0 


- 1 


2-1 

7 


(B2) 


If an asterisk is used to signify the ratio of the indicated parameter 
to the value of that parameter when the ratio of specific heats 7 is 
taken as constant, then the following relations are obtained: 
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/O 

* 

1 - 

(vrl 

5/2 !♦(>/ 

T p,») 

r<v 

1 v j 

(Vo)* + 

(p( T p,«) 


(v p 0 r 

(>V P °)* /(V T o)* 


(B5) 


(B 6) 


where (B is the Sutherland constant and Tp^oo the static tempera- 
ture evaluated for constant 7. The much used coefficient in the 
linear formula for viscosity becomes 


C 


* 



(B7) 


Equation (BJ) is valid for all wall temperatures. The hypersonic inter- 
action parameter X then becomes (for all wall temperatures) 



By using 5>500° R for the value of the molecular vibrational constant 0 
and y = 7/5, figure 11 was constructed from these equations . For the 
range of stagnation temperature shown (1,000° to 5>000° R), the static 
temperature must be less than about ^>00° R for excellent accuracy. The 
Mach numbers at the top of the figure indicate the approximate limita- 
tions at integral temperatures imposed by the assumption of low static 
temperature . 

Now, if the ratio of specific heats y is assumed to be a constant, 
a simple solution for imperfect-gas effects may be obtained. The equa- 
tion of state for a gas at moderate pressures, when only the first two 
terms in a series representation are used, is 
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p = p(KT + Bp) 


where B Is known as "the second virial coefficient in the expansion 
and R is the gas constant . The differential equation for an adiabatic 
process, to the same approximation, is (p. 180 of ref. 25 ) 

7-1 dp = dTA p dB\ (B9) 

7 P T\ R dT / 

If Berthe lot 1 s value for B (p. 119 of ref. 25) is used: 

q RT C ( 6T 0 2 \ 

B = 55 571 l "F7 


and 7 is taken as constant, the exact solution of equation (B9) is 


7-1 


7-1 

*)’ 


/ T„\5 P 


0 


1 + SL W 

32 


T J p c 


3 - 27 + 2l/TcA 
' 32\T ro ) P c 


(BIO) 


or, since the deviations from the perfect gas are assumed to be small, 


T 

00 




gl 7-1 

32 7 



3 - 27 2 

7 32 P c 
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APPENDIX C 


SKIN FRICTION OF A FLAT-PLATE WING WITH DELTA PLAN FORM 


INCLUDING BOUNDARY-LAYER DISPLACEMENT EFFECT 


In this appendix a solution is derived for the skin friction with 
boundary-layer displacement effect on the flat-plate wing shown in the 
following sketch: 


The assumption is made that each filament in the stream and boundary 
layer over the wing remains unaffected by adjacent filaments, and thus 
that any section of this wing may be treated as a section of an unswept 
flat plate with the local chord the same. The average skin-friction 
coefficient for the delta-plan-form flat plate is, then, 


Flow 




or, in terms of the hypersonic interaction parameter. 



(Cl) 


where 



and 



turn. R 


l c 
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and R are the Reynolds numbers based on local chord and root chord, 
c r 

respectively, and C is the coefficient in the linear formula for 
viscosity. 


If Cp is the average skin-friction coefficient obtained from 

classical flat-plate boundary-layer theory (in effect, the Blasius 
solution) , in order to simplify the solution the following relations 
are assumed: 




(C2) 


= 1 + a'X 



(C5) 


The value Xq is the matching point between equations (C2) and ( C 5) where 
it can be shown that 



(C4) 


When X ^ X Q , equation (C2) can be substituted into equation (Cl) to 
obtain 


C F,A 




From appendix D (eq. (Dl8)): 


C F, A 



(C5) 


(c6) 




where 
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Then equation (C5) is 


= 3X C > 


■T,A 


X„ X 


dX 

iijz 


(C7) 


vhich, upon integration, becomes 


C F,A _ 6 b ,TT l/ 2 

r 5 ° 

C F,A 


(** * X» - £) < c8 > 


When X = Xq, equations (C2) and (C3) can be substituted into equation (Cl) 
to obtain 


4/ p°° -9/2_ pX 0 -5 _ 

C F,A = 4Xc r ( b 'J^ C F * dx + J_ ~ 


Xo _ __4 J 


p X dX + /_ C F X dX + a’ /_ C F X dX 


which is combined with equation (C6) and the result is integrated to 
give 

|M=l + |a^ r+ X 0r 5 /|-^-|^l.-l 5 \ (xc r S><o) (09) 

c F,a V 3 x 0 J/ V V / v 


or, by substituting equation (C4) into equation (C9), 

-fcArl! 

10 


-F,A 


2 ^ ,Xc r 


(a% r S X ) 


(CIO) 


Equation (C8) may also be written in terms of a 1 as 


-F,A _ 


C F,A 


,v W 2 


-¥K) 


1 Xc r = i^ (cn) 
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The value of a T over a range of wall temperature for both air and 
helium may be obtained from an integration of the results given in refer- 
ence 22. These values are given in the following table, where T w is 
an arbitrary but constant wall temperature and Tq is the stagnation 
temperature since the analysis is based on a Prandtl number of 1, though 
as an approximation Tq can probably be taken as the recovery tempera- 
ture for Prandtl number ^ 1. 


T 

7 = : 

1.4 

7 = 

I .67 

x w 

T o 

a T 

*0 

a ' 

*0 

0 

0.0983 

10.17 

0.1864 

5.36 

.2 

.1769 

5.65 

• 3531 

2.83 

.6 

.3846 

2.60 

.8206 

1.22 

1.0 

.6831 

1.46 

1.524 

.656 

2.0 

a 1.78 

a .562 

4.200 

.238 


^Prom an extrapolation of the values given in 
reference 22 . 


A plot of the values of a 1 presented above is given in figure 16. 

In the present case, for use with the delta-wing results presented 
in figure 12 , the wings were assumed to be insulated; thus a* = 0.68 
for air and Xq = 1.46. The maximum value of X Cr for the delta wings 

used in the present experiments was about 0.35> so that equation (CIO) 
with the third term neglected was used, or 


'F,A 


= C 


F,A 


1 + 1.02X n 


(C12) 


The value of ^^v/C was taken as 282 and as 1.64. The solu- 

tion based on consideration of a reverse delta wing was assumed to hold 
for a delta wing flying point forward. 
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APPENDIX D 

SKIN FRICTION OF A FLAT-PLATE DELTA WING WITH 
BOUNDARY-LAYER TRANSITION 


This appendix is concerned with an evaluation of the skin friction 
on a flat-plate delta wing where a portion of the wing is covered by a 
turbulent boundary layer. 


Effective Chord 


The average skin-friction coefficient for a longitudinal element of 


a flat plate may be written Cp <x 



7-. If a triangular flat plate .is 

cVn 

being considered, with dimensions as 
shown in the sketch, the effective 
chord for purposes of determining 
the average skin friction for the 
entire plate is defined as 


5 - (-2- f ” ^ is 

\ bc rJo ) 


-n 


(Dl) 


Equation (Dl) may be nondimens ionali zed by letting s' = s/b and 
c* = c/c 


f c r , and the result is 


c 

c^ 


1 — 

> r (c') n ds 

Jo 


-n 


(D2) 


Since c f = 1 - s ! equation (D2) may be readily integrated to obtain 


c _ 


2n - 1 
2n 


n 


(D3) 
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For various values of n, equation (D5) gives results as follows: 


n 

c 

c r 

2 

O.5625 

5 

.5785 

4 

.5862 

5 

.5905 

6 

•5955 

7 

•5955 


Consider now a triangular plate with a cutout, as shown in the 
following sketch: 



The effective chord in this case (where the skin-friction coefficient is 
based on the area of the wing without the cutout), when n = 2, is defined 
as 


c 





or, in nondimens ional form. 


c 



1 

4 


•Vb 


(c 


Joels’ 



(D4) 
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Substituting the variation of c with s into equation (D4) yields 

-]l/2 

ds' + 

-2 


-L = I, 

/ 

n ho/b 

\(2o . a \ B . + a" 

c r 4 

Jo 

b o\ c r ) 


_ “I 


1/2 

/ 

/ bC ° . (1 - S') 

ds 1 ? 


Cp(h - b Q ) 



Equation (D5 ) , upon integration, becomes 
c _ 9 b o 1 


l6 


b c o 


- a 




- a 


i 1 - t )(|) 


1/2 


(D5) 


(d6) 


Equation (D6) is indeterminate when cp = e; however, it is a simple matter 
to reintegrate equation (D5) with this restriction, and the result is 


jc_ = 1 

c r 4 




(D7) 


Laminar Skin Friction 

For the flat plates considered in the preceding section an effective 
chord was defined in such a way that the average skin-friction coefficient 
for the entire plate could be written as 


Cp - 


A 




When n = 2, substitution of equation (D7) gives 


Cp = 


2A 

^ 1/2 



(D8) 
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The ratio of the skin friction for the plate with the triangular cutout 
(where the edges of the cutout are parallel to the leading edges) to 
the skin friction of the complete triangular plate is as follows (from 
eq. (D8) ) : 



(D9) 


Skin Friction Due to Turbulent Boundary Layer 


The cutout line is now assumed to represent instead the line along 
which instantaneous boundary-layer transition occurs. It is further 
assumed that transition occurs at a constant distance from the leading 
edge (cp = e) . 


As a basis for estimating the skin friction in the triangular 
region with turbulent flow, the momentum loss for the laminar boundary 
layer at the transition line is equated to the momentum loss for the 
turbulent boundary layer at the transition line. For a flat plate this 
is equivalent to a simultaneous solution of the equations for the momentum 
thickness. For the laminar boundary layer the momentum thickness may be 
represented by 



(DIO) 


and for the turbulent boundary layer. 



(Dll) 


where Axij is the streamwise distance from the calculated start of the 
turbulent boundary layer at zero thickness to the point where transition 
occurs. The form of equation (Dll) was actually obtained from considera- 
tion of the local skin-friction coefficient for turbulent flow in the 
form 


C f = C 


^ \^ n 
Poo u oo x j 


(D12) 
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Equating (DIO) and (Dll) gives 



m 


(D13) 


where ARp is the Reynolds number from the calculated start of the 
turbulent boundary layer at zero thickness to the point where transition 
occurs, and R c ^ is the Reynolds number based on root-chord length. 

Assuming that the average turbulent skin friction may be represented 
by an equation of the form 


where q is the value of n obtained for the turbulent boundary layer. 
Then if the skin-friction coefficient is based on the total plan-form 
area of the triangular flat plate, substitution of equation (D 13 ) into 
equation (Dl4) yields 


The skin-friction coefficient that can be attributed to the area with 
turbulent boundary layer, including the imaginary area defined by ARrp, 
is, when based on the total plan-form area of the triangular flat plate. 



(D14) 


AC F,turb = D ( 1 “ a ) - 


A m - 1 
C m 



(D15) 


turb 


D 


(1 - a) 




where t is a measure of the average chord in the region bounded by the 
lines for boundary- layer transition. When equation (D13) is substituted 
into this equation the result is 
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1+0 


c F,turb 


D(1 - a) 


r 1 
2 

r 

A m - li 
LC m 1 

Kr ) 1/2 

m 

m-1 

Rcr 



+ £ ( 1 - a) Rc 


l/l 


(D1 6) 


The average skin-friction coefficient for the triangular flat plate which 
is alined with the stream, with the assumptions of boundary-layer transi- 
tion that have been made, is 


Cp 



°F, turb " AC F,turb 


(D17) 


Laminar Turbulent 

where F c /F t is given by equation (D9), C Fjturb by equation (Dl6), 
and AC F>burb by equation (D15) . 


Evaluation of Constants 

A proper evaluation of the constants involved is the crux in deter- 
mining accurate values for the average skin-friction coefficient by this 
method. 

The constant A in equation (D17) may be evaluated from any of the 
available theories for the laminar boundary layer. For this report, the 
method given in reference 2 for hypersonic flow over sui insulated plate 
with a Prandtl number of 1 was arbitrarily chosen. Equation (15) of 
reference 2 was, however, modified by the result of equation (D3) given 
here so that 


Cp 


1 + 

3 



k 

3 




(D18) 
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1/2 
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and 9/6 is the ratio of the momentum thickness of the boundary layer 


to the total thickness of the boundary layer (linear velocity profile) . 
Values of (Q/b)^ as a function of ^ are given in reference 2 or 
they can be obtained from the approximation 


With the free-stream temperature given in °R, the Sutherland con- 

stant p has been taken as 177° F for hypersonic Mach numbers under 
ordinary wind-tunnel conditions. For lower supersonic Mach numbers and 
for flight conditions, the Sutherland constant may be taken as I 98 0 F. 

For the evaluation of the turbulent part of the skin friction, 
values must be determined for C, D, m, and q. To this end Van Driest 's 
analysis for the turbulent boundary layer on a flat plate (ref. 26 ) was 
arbitrarily chosen; of course, there are other turbulent boundary-layer 
analyses which could as easily be utilized. The values of m and C 
were obtained by fitting equation (D12) in the range of Reynolds numbers 
desired to the values of local skin-friction coefficient obtained from 
equation ( 66 ) of reference 26. This equation is based on the power law 
for viscosity; however, for high temperatures where the variation of 
dynamic viscosity with temperature departs markedly from a power-law 
variation, it is a relatively simple matter to correct the exponent in 
the power law according to Sutherland's viscosity formula to obtain 


where oj is the exponent in the power law and T^ is the wall tempera- 
ture. If the plate is insulated and the Prandtl number is 1, equa- 
tion (DI 9 ) may be written as 


(§). 


O. 3 O 683 
2 + k 



(D19) 



log 10 ( 1 + k) 


(D20) 
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The values of D and q introduced in equation (Dl4) cannot, in 
this case, necessarily he found to the required accuracy from an inte- 
gration of equation (D12) . For the present computations, equation ( 79 ) 
(or eq. ( 77 )) of reference 26 was used to obtain values to which equa- 
tion (Dl4) was fitted over the desired range of Reynolds numbers. 

A value must now be assigned to |. This constant is actually the 
ratio of the effective chord in the triangular region with turbulent 
boundary- layer flow to (l - a)c r . When a = 0, ARip will also be zero 
and the value of | is that given by equation (D 3 ) for the value of 
n (=q) involved. A value of O .55 was used for t in the computations 
for this report, since the turbulent boundary layer assumed does not 
start with zero thickness. 


The remaining constants used were as follows: A = 1.22, C = 0.125, 


D = 0.243, m = 3.27, q. = 3-03, 
were assumed. The product aR Cr 
designated as Rip in figure 15- 


o> = 0 . 86 , and values of a and R c 
(the transition Reynolds number) is 


«r 
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Figure 1.- Pressure distributions in the three-dimensional nozzle. 
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Figure 2.- Square -plan -form wings. 
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Semi apex 
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in. 
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30 
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20 
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13 
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Figure 3 . - Delta wings . 
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+ Pressure orifices 



Pressure-Orifice Locations 
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Figure 4.- Flat-plate pressure and temperature model 
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Figure 5*- Theories for predicting the induced-pressure effects on an 

insulated flat plate, y = 7 / 5 . 
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Figure 6.- Theories for predicting the induced-pressure effects on local 
and average skin-friction coefficients for an insulated flat plate . 
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Figure 7 . - Comparison of induced pressures on a flat plate as given by theory and by 

experiments in air at M = 9.6. 
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Figure 8.- Measured wall temperatures on pressure model and prediction 
of Chapman -Rubesin theory for the temperature effect on the displace- 
ment thickness. M = 9-6. 
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Figure 9 .- Lateral distribution of pressure at two stations on flat 

plate. M = 9 . 6 . 
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Figure 10.- Induced pressures on insulated flat plate. Data for 
M = 9-6 with heat transfer but corrected theoretically to 
insulated-plate case. 
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Figure 11.- Effect of caloric imperfections (variable y) on 
Reynolds number , viscosity coefficient, and the viscous - 
interaction parameter. 
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Figure 12.- Experimental and theoretical drag coefficients of square - 

and delta -plan -form wings . 
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Figure 13.- Surface oil-flow studies of the square -plan -form wings. 

M = 6 . 9 . 
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(c) Wing 4-B; € = 8 °; R = 3.4 X 10 6 . L -87898 

Figure 14.- Surface oil-flow studies of two of the 2^ - percent -thick 

delta wings. M = 6 . 9 . Lines have been added to indicate edges of 
wings, which were not visible in photographs. 
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Figure 15.- Comparison of laminar and turbulent skin-friction coefficients with experiment. 
Experimental values are total drag from which an estimated leading -edge drag and pressure 
drag have been subtracted. M = 6.8. Open symbols denote data obtained in steel nozzle j 
solid symbols denote data obtained in Invar nozzle. 
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Figure l6.- Values of the coefficient a' in the equation for skin 
friction from viscous -inter act ion theory. See appendix C. 
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